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Abstract 

The Earth, like other planets with a relatively thick atmosphere, 
is not locally in radiative equilibrium and the transport of energy by 
the geophysical fluids (atmosphere and ocean) plays a fundamental 
role in determining its climate. Using simple energy-balance models, 
it was suggested a few decades ago that the meridional energy fluxes 
might follow a thermodynamic Maximum Entropy Production (MEP) 
principle. In the present study, we assess the MEP hypothesis in the 
framework of a minimal climate model based solely on a robust ra- 
diative scheme and the MEP principle, with no extra assumptions. 
Specifically, we show that by choosing an adequate radiative exchange 
formulation, the Net Exchange Formulation, a rigorous derivation of 
all the physical parameters can be performed. The MEP principle is 
also extended to surface energy fluxes, in addition to meridional energy 
fluxes. The climate model presented here is extremely fast, needs very 
little empirical data and does not rely on ad hoc parameterizations. 
We investigate its range of validity by comparing its performances 
for pre-industrial climate and Last Glacial Maximum climate with 
corresponding simulations with the IPSL coupled atmosphere-ocean 
General Circulation Model IPSL_ CM4, finding reasonable agreement. 
Beyond the practical interest of this result for climate modelling, it 
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supports the idea that, to a certain extent, chmate can be character- 
ized with macroscale features with no need to compute the underlying 
microscale dynamics. 

1 Introduction 

The Earth receives a certain amount of energy from the Sun, in the form of 
visible light, which it has to radiate back to space, in the form of infrared 
light, to maintain a steady state. Most Earth System processes, including 
weather and climate, can be regarded as little more than steps in this process 
of energy conversion from one form to the other, going through various other 
forms of energy (potential energy, kinetic energy, heat...). Although in each of 
these steps the overall quantity of energy has to be conserved, its quality can 
vary |^ HQ]. The field of physics that deals with such energy conversions 
in different forms and the corresponding evolution of its quality is called 
thermodynamics, and the quality of energy is measured by a state function 
called entropy. The constraint of conservation of the exchanged quantity of 
energy is well-known and commonly used in climate sciences. Less attention 
is paid to the additional information one can gain by monitoringcreation and 
exchanges of entropy, either as a diagnostic tool or as a prediction principle 
relying on an extremum property. This information can be discarded with 
little consequences for some processes for which physical laws naturally allow 
for a macroscopic description (and for the subsequent computation of the 
exchanged energy due to the process), but in some other cases it might bring 
in a fundamental constraint that should not be ignored. 

A typical example is atmospheric heat transport: considering the obser- 
vational fact that the net top-of-atmosphere radiative budget is positive in 
the Tropics and negative near the Poles, one has to invoke a poleward at- 
mospheric transport to ensure global energy balance in the steady state (see 
e.g. [50j for a quantitative discussion based on observations and [28j for a 
modelling counterpart). An important difficulty comes from the turbulent 
nature of the laws of atmospheric motion and the lack of a satisfactory theory 
of turbulence. We are thus forced to integrate numerically the equations of 
motion at a very high computational cost. On the other hand, the laws of 
radiation allow for a macroscopic description more easily. To cope with this 
difficulty, the options up to now have been either to deal with the overdetailed 
(from the minimalist standpoint we are considering) microscopic equations 
of motion of the atmosphere, or to use empirical parameterizations. 

An alternative was suggested by Paltridge [35l l36] (see |33) for a nice 
reformulation). He was looking for a large-scale description of heat transport 
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and tested a variational principle stating that atmospheric heat transport 
adjusts so as to maximize the production of entropy. Several climate studies 
based on similar models have been published since this seminal work: [T3| 
[TT| inU |25] for zonally-averaged climate, [Ml H3] for vertical convection. The 
Maximum Entropy Production can be interpreted as a maximum energy 
transport efficiency requirement. The general idea was first expressed by 
Lorenz |2l]. See for instance |10t |23| |52| [26l [27] and references therein for 
discussions in terms of available potential energy and Carnot engines, and 

IBS] for the competing effect of water vapour removal by the atmosphere. 

The purpose of this paper is not to discuss the justification of the Max- 
imum Entropy Production Principle. For a rather comprehensive review of 
different approaches and applications, the reader is refered to [31j and |20] . 
In this paper, we present a new formulation of Paltridge's model: we have 
developed a different approach for the treatment of radiation, which we be- 
lieve to be more rigorous and physically sound. In particular we provide 
a derivation of the analytical expression of the radiative coefficients which 
are introduced as ad hoc in Paltridge's model. We also avoid the question- 
able maximum convection hypothesis, critical in Paltridge's model [33| HI], 
by applying the MEP principle on the vertical dimension as well, and dis- 
card Paltridge's cloud cover variable and oceanic transport for the sake of 
cogency. Formulated this way, the model is as free of tunable coefficient, 
empirical parameterization and spurious assumption as possible, therefore 
constituting a clean basis to evaluate the validity of the MEP hypothesis as 
applied to climate. The only information needed a priori to compute the co- 
efficients are the surface albedo, the vertically-integrated water vapour and 
ozone density and the CO2 concentration. In particular, we are now able to 
run a sensitivity experiment with respect to the surface albedo parameter. 
Here, we investigate the effect of the surface albedo distribution correspond- 
ing to the presence of large ice-sheets in the Northern Hemisphere during the 
Last Glacial Maximum. The results for pre-industrial and LGM climates are 
compared with standard simulations with the IPSL_CM4 atmosphere-ocean 
general circulation model using similar boundary conditions. 

2 Model description 

We use a purely thermodynamic two dimensional (horizontal) energy bal- 
ance model. Only energy transfers are represented. The exchanges of mass, 
momentum, angular momentum, etc are discarded here. Each grid cell con- 
sists of two layers: one atmospheric layer (at temperature T^) exchanging 
energy with the surface, space and adjacents cells in the form of longwave 
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and shortwave radiation, surface heat flux, and atmospheric heat transport, 
and one surface layer (at temperature Tg), only exchanging energy with the 
overlying atmospheric layer and with space, in the form of shortwave and 
longwave radiation and surface sensible heat (see figure [T]). We note Aij the 
area of cell 1 < « < Niat, 1 < J < iV/on- All the simulations presented 

in this paper are carried out on a regular grid with Niat = 72 and Nion = 96 
for easier comparison with the GCM results using the same grid. 

2.1 Net-exchange formulation for longwave radiative trans- 
fer 

Radiative exchanges are described in the Net Exchange Formulation (NEF) 
formalism - originally introduced by |14j (see [3], [7] for details). In the NEF, 
radiative fluxes are replaced by energy exchange rates between two layers 
as basic objects describing the radiative transfer. Instead of solving the 
radiative transfer equation, we compute the net energy exchanged per time 
unit between layers i and j (resp. layer i and boundary surface j or boundary 
surfaces i and j) at frequency u: 

III dPdP.d^ 

Jt,, JSj Jr., (1) 

where Ej represents the volume of layer i (resp. the area of boundary 
surface i), Fjj the set of optical paths from one point of layer i to one point of 
layer j, aj,(Pj,7) is 1 if i is a boundary surface and the absorption coefficient 

k^{Pi, 7) otherwise, = exp ^— k^{s)ds^ is the transmitivity of the path 

7 and finally By{T) is the Planck function. Depending on the nature of Ej 
(surface or volume), dPi is either an infinitesimal surface or volume element. 
Clearly r^, = -V.i- 

The NEF was originally used in engineering rather than atmospheric sci- 
ences. In the few existing cases of atmospheric transfer studies, the main 
point is that NEF allows for a fine analysis of the greenhouse effect and how 
it evolves when one parameter changes: it is clear from the {"ipij) matrix which 
exchanges prevail in the atmosphere (generally, these are the exchanges be- 
tween neighbouring layers, between the surface and any layer and the cooling 
to space of any layer and the surface). This description contains more infor- 
mation on the radiative transfer than the traditional fiux approach insofar 
as one not only knows how much energy is radiated through a surface but 
also where the radiation originates from (computationally, for a n layer at- 
mosphere, NEF requires 0{v?) numbers while the fiux approach only needs 



4 



0{n)). Other major advantages of this formulation are the handiness it pro- 
vides in couphng radiative exchanges to other modes of energy transfer, and 
the fact that it automatically enforces physical laws such as energy conser- 
vation or the direction of energy exchanges depending on the temperature 
gradient, whatever numerical approximations are made. Clearly these points 
make it suitable for simple thermodynamic models such as ours. 

In the very simple case of our two level cell, we only have to compute the 
infrared exchange between the surface and the atmosphere ^E'^f^, the surface 
and space "^^f and the atmosphere and space ^E^^f . Applying the diffusive 
approximation (/i = 1/1.66 [4J, justified in |22)). angular and vertical in- 
tegration [H being the height of the atmosphere) yields, for instance, for 
•qj" ■ 

0-9' 

= (l - e-fJ"''^^^^'^^ {B,(Tg) - B,{T,)) (2) 

From now on, let us consider the exchanged energy per surface area, 
keeping the same symbols \E'y. Integrating over the whole spectrum, the 
equations for the exchanged energy can be recast into the following form: 

^af = UTg)aT'^-UT,)aT', (3) 
^if = t^iT^)cyn (4) 
K = (5) 

1 — / — T^du ] represents the emissivity of the at- 
Jo crT4 J 

mosphere and t^, = exp ^— ^ f^^ k^{z)dz^ . Of course the transmission func- 
tions of the atmosphere depend on the vertical profiles of absorbing gases, 
pressure and temperature. In order to preserve the simplicity of the analyt- 
ical expressions and to allow for comparison with Paltridge, we impose here 
the heavy constraint that the coefficients ta are constant in each cell. Given 
a characteristic temperature for the local atmosphere and a humidity profile, 
we can evaluate these coefficients using Goody's statistical model [12j and 
the spectral data from |45| . The temperature and humidity profiles used here 
come from a linear interpolation of McClatchey's reference atmosphere pro- 
files [32]. Hence we prescribe here the ta coefficients as functions of latitude 
only. In fact we do not need the full profiles here but only the vertically inte- 
grated values. Sensitivity tests were carried out to assess the dependence of 
the ta coefficient on the water vapour amount and characteristic temperature 



and its impact on the MEP climate (see section 3.3). 
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2.2 Shortwave radiation parameterization 

At latitude 0, the annual mean insolation is given by 

S 



Wa 



"27r2x/r^ 

/•27r 

/ a/1 — (sin (j) cos e — cos sin e sin u^du 
Jo 



(6) 



where is the solar constant, e the eccentricity of the Earth's orbit and 
e the obliquity. The incoming top-of-atmosphere solar radiation in each cell 
is given by the latitudinal mean of Wa over the cell, which we write ^ x S*, 
^ thus representing the annual mean of the cosine of the zenith angle of the 
Sun for the latitude zone (weighed by a distance factor) and is therefore a 
function of the orbital parameters. In all the simulations considered here, 
we use the present-day values for the orbital parameters: e ~ 0.0167 and 
e ~ 23.4 ° . The Last Glacial Maximum orbital parameters do not differ 
much from those values. 

Solar radiation exchange rates are parameterized as in [2T]. With their 
notations, let us define 

Sa{Rg) = 0.353 H = (7) 

1 — R*Rg 

Sa = A^.,{Mu) (8) 
< = A^, (^(^M+^^u^- A^,{Mu) (9) 

where Rg is the surface albedo, u (resp. Uo-j) the column-integrated water 
vapour (resp. ozone) concentration. M is the magnification factor accounting 
for the slant path, Ay^^ (resp. Aoz) are numeric functions j3orresponding 
to water- vapour (resp. ozone) absorption, and i?r(0 R* account for 
Rayleigh scattering in the atmosphere (see [21j or [46j for details). Incoming 
solar radiation at surface and solar radiation absorbed in the atmosphere are 
then expressed in the simple form: 

= (SaiRg) - S,){1 ~ Rg)^S (10) 

v^ff = {s^ + RgsD^S (11) 

The coefficient Sa can be interpreted as the proportion of solar radia- 
tion directly absorbed in the atmosphere while s* accounts for absorption 
after reflection at surface and Sa represents the proportion of solar radiation 
absorbed by the ground. 
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To compute the value of the three coefficients Sq, Sa, s*, we use humidity 
and ozone profiles from McClatchey |32| again (as for the longwave coeffi- 
cients, this is not a critical parameter here), with pressure scaling for water 
vapour. Surface albedo is prescribed from IPSL simulations (see section 2.4). 



2.3 Entropy Production Maximization 

For each grid box, we can compute the parameters (^, t^, Sq, s*, s^, i?^) as 
explained above and radiative energy exchange rates are given by equations 



|3]to 11 as functions of these parameters and temperatures Tq and Tg. 

Since we are looking for steady states of the system, we impose the fol- 
lowing energy balance constraints in each box: 

^r+^ff-^^f-^- +c = (12) 

where q is the surface heat flux and ( the horizontal convergence (see 
figure [T]). Formally these equations can be solved to express Ta and Tg as 
functions of the two unknown variables, q and (. 

Similarly to the energy budget, we can write an entropy budget, sepa- 
rating the rate of change of the entropy of a parcel of the system into two 
contributions: 

dS dSe 

— = —^ + a 14 
dt dt ^ ' 

where Se stands for the entropy exchanged with the surroundings per unit 
time, while a is the entropy created per unit time. The second law of ther- 
modynamics states that a > 0. In the steady-state, we have S" = 0. 

Assuming local thermodynamic equilibrium, the variation of the entropy 
for a parcel of the system exchanging an amount of heat 5Q is given by 
dS = 6Q/T. Hence, for a single cell of the model, the entropy production rate 
(per unit surface area) associated to the horizontal convergence is given by — 



and the entropy production associated to the surface heat flux is q 
We thus consider the total material entropy production rate: 



Ta 

J 1_ 

Ta T„ 



(^M{{qijXij]) 



Nlat Ni^n 

EE 



i=l j=l \ o^'^J 9,1-J -^a,!] 



Let us search for the energy fluxes distribution maximizing this function 
subject to the following global constraint: the total heat transport divergence 
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over the globe must vanish, i. e. 

Nlat Nl^n 

^J]A,,O, = (16) 
i=i j=i 

We enforce this physical condition by introducing a Lagrange muhipher /3 in 
the previous equation. 



<^M{{qij,Cij},^) 

Nlat Nl^n 

i=l 3 = 1 



T, 



a,ij 



T, 



+ 



/3C 



(17) 



and we call Maximum Entropy Production (MEP) state the solution 
{^a,ij.^g,y.9*i.Q}i<i<^^^^j.<^^ of the system (|l2|)-(|l3l maximizing the 



entropy production rate (15) with the imposed constraint (16). In practice 



solving the system reduces to finding the maximum of a function of two vari- 
ables in each box. One could equivalently solve the optimization problem in 



terms of the temperature variables, replacing q and C, from equations 12 13 
into (TM{Ta,Tg). 

In the MEP procedure, it is only the entropy production rate due to tur- 
bulent processes which is maximized, motivated by the supposed connection 
of the MEP hypothesis with principles of maximum statistical entropy. The 
entropy production rate due to radiative processes, however, can provide 
interesting information as a diagnostic tool (see section 3.5). 



2.4 IPSL simulations description and albedo forcing 

The IPSL model is a fully coupled Ocean-Atmosphere General Circulation 
Model (OAGCM). The version used here is IPSL_CM4 [3D] i.e. the version 
which has been used for the CMIP3/IPCC AR4 exercise. The model is 
composed of the LMDZ atmosphere model [16] run at resolution 96x72x19 
(longitude-latitude- vertical levels), the 0RCA2 ocean model [29j . the LIM2 
see-ice model [HIE], all coupled via the OASIS coupler [51| . 

The present work uses two simulations performed with this model : (1) 
a pre-industrial simulation, for which the model is forced with constant, 
pre-industrial greenhouse gas concentration {CO2 = 280 ppm, CH^^ = 760 
ppb, N2O = 270 ppb), insolation, coastlines, topography, land-ice extent 
and (2) a pseudo LGM run. Since the simple climate model described above 
does not account for topography, nor for CO2 variations, we have chosen 
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to use an OAGCM simulation in which the extent of the LGM ice-sheets is 
taken into account but not their height. The ice-sheet extent is prescribed 
from the Peltier ICE-5G reconstruction [42j. Greenhouse gases atmospheric 
concentrations are maintained at their pre-industrial level, as well as orbital 
parameters and topography. Note that none of the simulations uses an aerosol 
forcing, which explains that the simulated pre-industrial climate is slightly 
warmer than the standard IPCC simulation. 

For each of these simulations, the annual mean surface albedo is computed 
(figure [2]) and used as an input for the corresponding MEP simulation. 

Hence the boundary counditions for the IPSL_CM4 simulations corre- 
spond as much as possible to those for the MEP simulations. 



3 Results and discussion 
3.1 Present day climate 

Figure |3] shows the surface temperature at MEP state for present albedo 
forcing and a comparison with a control simulation from the IPSL model. 

The global mean temperature simulated by our simple MEP model is 
23 ° C, which is 7 ° C warmer than the IPSL model results. The reason for this 
strong bias is probably the lack of any clouds in the model. Another reason 
could be the crude representation of the vertical structure of the atmosphere 
in the model and the subsequent flaws in vertical radiative transfer calcu- 
lations. Comparison with the IPSL model reveals that some of the regions 
of largest discrepancies are regions of high elevation (Antarctica, Tibetan 
plateau, and to a smaller extent Greenland, cf Fig. Isj). Indeed, topogra- 



phy effects are not included in our model. We show in section 3^ that the 
combined effects of clouds and elevation account for the bulk of the global 
difference with the IPSL model. 

In spite of this global bias, considering the simplicity of the model and 
the absence of any degree of freedom to adjust the results, the model gives 
a strikingly accurate view of the global picture of the climate. Some impor- 
tant specific features, like the presence of the main deserts (Sahara, Arabian 
peninsula, Kalahari, Australia,..), are already visible taking into account only 
the albedo effect. Since the only non zonal parameter in the model is the 
albedo, and its value is roughly uniform over the oceans, we were expecting 
the surface temperature to be zonally homogeneous over the oceans, which 
is indeed the case. The positive anomaly of temperature over mid-latitudes 
regions (compared to the IPSL results. Fig. [sl right) is certainly imputable 



to the important cloud cover over those areas, lacking in the model. 
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Figure |4] provides the curves of zonally averaged atmospheric heat trans- 
port and net top of atmosphere radiative budget. The net top of atmosphere 
radiation budget switches from positive in the Tropics to negative near the 
Poles at roughly 40 ° S and 35 ° N, while the overall shape of atmospheric heat 
transport qualitatively fits the IPSL model results. It peaks approximatively 
at 35 ° S and 40 ° N with respective values of 3.8 PW and 4.4 PW. The loca- 
tion of these maximum values is consistent with |l8] but the absolute value 
is slightly underestimated here: gives the peak transport around 6 PW. 
The absence of any oceanic heat transport in our model prevents us to push 
the comparison forward. Comparison with the atmospheric transport from 
the IPSL simulation reveals that the MEP transport is indeed a little low, 
especially in the Southern Hemisphere. 

It might also explain why surface energy fluxes are underestimated here: 
since all the heat flowing laterally is assumed to come from the atmosphere, 
the need for surface heat flux is reduced. Indeed, in the MEP state, the global 
mean surface flux is only 55W.m~^ (figure [5|, compared to 97W.m~^ |49) . 
The maximum value of the surface energy flux is reached in the Tropics while 
it is negative in the polar regions. Comparison with the IPSL simulation (not 
shown) confirms that the largest discrepancy arises over the oceans, where 
the lack of oceanic transport in our MEP model is crucial. 

These results can also be interpreted in a different manner, considering 
that the surface energy flux computed by MEP is actually closer to the 
sensible heat component rather than the sum of the sensible and latent heat 
flux. Formally, the physical nature of the surface flux is not specified in the 
equations of the model. Yet, we believe that since it appears in the entropy 
production rate as the product between the flux and the gradient of the 
inverse of temperature, and given that the sensible heat flux also depends 
directly on the temperature gradient at the surface but not the latent heat 
flux, it makes more sense to compare it to the sensible heat flux (figure |5]). 
In that case, the global average surface flux is approximately twice the value 
given by [49j. The fact that the global mean average surface flux lies in 
between the value of the global mean sensible heat flux and the sum of the 
sensible and latent heat flux can be seen as a trade-off between the need 
for a large surface flux to ensure the energy balance and the need for a large 
enough temperature gradient from the entropy production rate point of view. 

As discussed in jl9], since a MEP state is supposed to be the most efficient 
heat conducting state given the constraints, one should expect more gentle 
equator-pole temperature gradients in MEP states than with more conven- 
tional models. Here, adjusting the overall surface temperature field by setting 
Tg = Tg + {{Tg)ipsL — {Tg)MEp) to eliminate the global bias, we observe that 
the Tropics are roughly at the same temperature in the IPSL model and 
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the adjusted MEP model (actually slightly colder in the MEP model) while 
the mid-latitudes are substantially warmer (around 5 degrees) in the MEP 
state. However the temperature anomaly is not symmetric for the poles: the 
South Pole is colder by around 10 ° C in the adjusted MEP model whereas the 
North Pole is warmer by around 5 ° C. Therefore the equator-pole gradient is 
indeed reduced in the MEP state for the Northern Hemisphere, but not for 
the Southern Hemisphere, as expected from the relatively low poleward heat 
transport in the Southern Hemisphere (figure ^ . 



3.2 A representation of clouds in the MEP climate 

As suggested in section |3.1[ the major part of the difference between the 
MEP pre-industrial climate and the corresponding IPSL simulation is likely 
to be due to the absence of clouds in the MEP model. Inclusion of fixed or 
variable clouds in MEP models is possible, but unavoidably requires param- 
eterization and extra assumptions, which is not the philosophy of this study. 
Nevertheless, to evaluate the effect of clouds on the MEP model presented 



here, it is possible to modify the radiation equations (lO)-(ll) to include an 
ad hoc cloud parameter 



= (SaiRg) - Sa){l - Rg){l - Rn)^S (18) 
= {Sa + RgS:){l-R^)^S (19) 

The parameter Rn is computed in each grid box from the cloud radiative 
forcing at the top of the atmosphere in the corresponding area in the IPSL 
model. On the global mean, the cloud radiative forcing at the top of the 
atmosphere in the IPSL model is roughly —20W.m~'^ (in accordance with 
|18)). which corresponds to a mean value of i?„ ^ 0.073, but the geographical 
distribution of the radiative forcing is very far from being uniform. 

For pre-industrial climate, this parameterization of clouds in the MEP 
model reduces the temperature by 4.7 ° C on global average. Thus a 2.5 ° C 
difference remains with the corresponding IPSL run. The temperature anomaly 
map (not shown) reveals that this difference is mainly due to the effect of 
topography (Antarctica and Tibet are in particular very warm in the MEP 
model), not accounted for in the MEP model. As an approximate correction, 
we apply a constant —7K.km~^ lapse rate, using the elevation data from the 
IPSL run output. Figure |6] shows the combined effect of clouds and topogra- 
phy on the MEP pre-industrial climate: the global mean difference with the 
IPSL pre-industrial run becomes 0.9 ° C, which is comparable to the uncer- 



tainty due to the water vapour content of the atmosphere (see section 3.3). 
Thus a crude parameterization of the main components lacking in the rigor- 
ous version of the MEP model suffices to recover the correct global picture of 
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present climate. Note in particular that the warm bands in the mid-latitudes 
in figure [3] have completely disappeared with the inclusion of clouds. 

3.3 Effect of the radiative parameters on the MEP state 

As mentioned above, the only parameters in the model are the radiative pa- 
rameters {ta, Sa, s*, So) and the surface albedo Rg (forced in each box with the 
surface albedo from the IPSL run). Computation of the radiative parameters 
only requires the value of the column-integrated water vapour density UH2O) 
and a characteristic temperature T for ta- 

In this study, we have fixed once and for all the value of the radiative 
parameters as a function of latitude only, using a linear interpolation of 
standard atmospheric profiles for UH2O and T. To assess the sensitivity of 
the model with respect to this choice, we have also computed the MEP state 
for pre-industrial surface albedo using a constant value for t^, independent 
of latitude. Three choices were made, corresponding to the three standard 
profiles: Sub-Arctic, Mid-Latitude and Tropical. Table [T] shows the global 
mean temperature resulting from these uniform choices. 

The values for the column-integrated water vapour density and charac- 
teristic temperature span a realistic interval: as a comparison, in the IPSL 
pre-industrial run, UH2O ranges from around 0.25g.c'm~'^ near the poles to 
around b.bg.cm^'^ at the equator. Thus we can estimate the sensitivity of 
the model to realistic choices in UH2O and T to be on the order of IK. 

3.4 Last Glacial Maximum 

Figure [7] shows the surface temperature difference between the Last Glacial 
Maximum and pre-industrial control run, at MEP and with the IPSL model. 
The global mean difference is ^ 2 K colder for the LGM at MEP, compared 
to ~ 2.5 K colder for the IPSL model. This remarkable overall agreement 
means that the sensitivity of the MEP model to surface albedo is comparable 
to that of the GCM. The major part of the cooling occurs over Canada and 
Northern Europe for both models, but it is much stronger in the MEP model: 
between 20 and 30 K for these areas, compared to ~ 15K for the IPSL model. 
The cooling is also overestimated in Patagonia - by a factor 2 approximatively 
- and Tibet - by a factor 3. These extreme values are compensated for by 
a small underestimate of the cooling almost everywhere else: around —0.5 
K for most of the oceans. South America, Africa, Australia, Antarctica and 
Southern and Eastern Asia vs — 1 to —2 K for the IPSL model. This milder 

cooling becomes even more dramatic over Russia 2 K for MEP vs —6 K 

for the IPSL model - and Greenland — 1 K vs —5 K. Surprisingly the MEP 
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model predicts an intense warming - ^ 8 K - in the Baffin Bay. We suspect 
this warming to be due to seasonahty effects of the albedo. 

Beyond the general relative agreement, the MEP model response to sur- 
face albedo is a little too strong over land (Canada, Northern Europe, Patag- 
onia and Tibet) but not over sea: the Barents sea is slightly cooler in the 
IPSL-LGM than in the MEP-LGM. 

Figure |8] shows the meridional heat transport for the Last Glacial Maxi- 
mum in MEP state. There is no major change of regime compared to present 
climate. Poleward transport is slightly more important (roughly 10% more 
energy at the maximum) during the LGM, as can be expected from colder 
high latitudes. This enhancement of the MEP-LGM poleward transport, 
larger than that obtained by the IPSL simulation, can probably be attributed 
to the exaggerated cooling of the Northern high latitudes in MEP-LGM. Dur- 
ing the Last Glacial Maximum, surface fluxes (not shown) become negative 
over ice-covered areas, the remaining of the globe showing little difference 
with present surface fluxes at MEP state. 

The effect of clouds for the LGM at MEP state, with the cloud parameter 
being calculated from local cloud radiative forcing in the IPSL LGM simula- 
tion, is less than for pre-industrial climate. On the global mean, clouds cool 
the LGM climate by 4 ° C in the MEP model. As a result, the sensitivity of 
the MEP model with clouds is reduced to — 1.3 ° C, leading to a sensitivity 
difference with the IPSL model of 1.2 ° C, as compared to 0.5 ° C for the 
model without clouds. Topography has no influence on the LGM-PI surface 
temperature difference as it is the same for both periods in both models. 

Since the radiative coefficients are set once and for all in the MEP model, 
there is no representation of the water-vapour feedback, contrary to the 
IPSL model. In the IPSL model, the globally averaged water-vapour den- 
sity changes from 2.86g.cm~'^ for pre-industrial climate to 2.50g.cm~'^ at the 



LGM, corresponding to a variation in of around 0.003. From section 3.3 
one can estimate that the effect of water-vapour feedback on the MEP-LGM 
climate must be on the order of a few tenth of degrees. 

3.5 Entropy budgets 

We can decompose the total entropy production rate of section 2.3 into con- 
tributions from various processes as follows: 

^ ntm I sfc I atm—sfc . /'on^ 

^ = ^SW + ^SW + ^LW + (20) 

where cr^iy , ct^vk, cr'^w~'^^'^ entropy production rates due to, respec- 

tively, absorption of solar radiation in the atmosphere, absorption of solar 
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radiation at the surface and surface-atmosphere longwave interaction: 



atm—sfc 
LW 



a 



.atm 
SW 



.sfc 
SW 




(21) 



(22) 



(23) 



(24) 



The entropy production rate of a thermodynamic process can be inter- 
preted as a measure of its irreversibihty. Therefore, the inspection of the 
relative values of the different terms above, gives valuable insight with re- 
spect to the role played by the different processes in driving the Earth out- 
of-equilibrium. 

The expressions for the contributions to entropy production involving 
interaction between radiation and matter are only approximate, as noticed 
by [5], but we shall not delve further into this here, both for simplicity and 
to compare with published results. The reader is referred to [47j or |53| . 
for example, for computations taking into account the correct entropy of 
the radiation field. Nevertheless, the nonlocality of the NEF is already an 
improvement, in principle, as compared to the usual flux formulation, which 
is a local expression of radiative transfer [6j. 

Table |2] shows the entropy budget of the Earth under pre-industrial and 
Last Glacial Maximum conditions at the MEP state. These values are com- 
pared with estimations for present climate from observations [H], from the 
intermediate complexity model PLASIM |10) and from the Hadley Centre 
GCM |37j. In both cases, the entropy production due to absorption of so- 
lar radiation by the atmosphere is a little lower than other estimates. More 
importantly, the entropy production due to absorption of solar radiation at 
the surface is significantly higher than other estimates: this is due to the 
high solar absorption at the surface (in the absence of clouds) which is not 
compensated for by a warmer surface. The entropy production due to long- 
wave interaction between the surface and the atmosphere (the greenhouse 
effect) is comparable to the values for PLASIM and HadCM3, and the ma- 
terial entropy production is considerably lower than other estimates. This 
last point is explained by the fact that latent heat fluxes are not explicitly 
represented in the MEP model, while they account for the major part of 
the material entropy production in all the other studies. Note that the total 
entropy production at MEP state is slightly lower for the LGM as compared 
to the pre-industrial, in spite of a higher material entropy production. 
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4 Conclusion 



The aim of this study was to provide a physically reliable, more consistent and 
more easily generalized version of the Paltridge model. The radiative part 
has been entirely reformulated using the Net-Exchange Formulation in which 
the physical meaning of all the coefficients is clearly identified. Therefore the 
calculation of these coefficients, to the desired degree of approximation, only 
involves standard radiative data. Besides, all the energy fluxes other than 
radiation (horizontal convergence and surface fluxes) are subjected to the 
same principle: Maximum Entropy Production. 

We have shown here that a thermodynamic model based only on a simple 
but robust radiative scheme and the principle of maximum entropy produc- 
tion yields results comparable to those of the IPSL model, for surface temper- 
ature, for both the pre-industrial and the Last Glacial Maximum climates, 
with a negligible computational cost. Obviously, the description adopted in 
the MEP model discards an important number of climate features. The in- 
clusion of some key processes, like the seasonal or the hydrological cycles, in 
this thermodynamic framework would be of great interest. But the model in 
its present form already serves the purpose of providing a compromise be- 
tween a minimal model with very little a priori data and information needed, 
and acceptable efficiency. 

In particular, one of the strongest points in the model is certainly that 
it does not include any adjustable parameter. This ability to get rid of the 
usual, varyingly important, parameter calibration makes it a good candidate 
to investigate climates where little is known or where some phenomena are 
likely to be different from the usual parameterization validity range, like 
for instance climates of other planets, inside or outside the solar system, or 
paleoclimates. 

Further developments also include coupling the model with a state-of- 
the-art radiative code to come up with a full three dimensional model, which 
would allow for more realistic paleoclimate simulations. 

The basic idea of variational thermodynamic principles such as the prin- 
ciple of Maximum Entropy Production, which still remains to be proved, is 
that microscopic details are irrelevant, to a certain extent and depending on 
the question we are trying to answer, to macroscopic behaviour. The gen- 
eral process of not making any unjustified assumption on those microscopic 
details finds a nice formulation at equilibrium with the principle of Maxi- 
mum Entropy and the usual apparatus of equilibrium thermodynamics can 
be derived easily from the application of this principle to the classical ensem- 
bles of statistical mechanics [17j. But in spite of the conceptual similarity, 
the link between Maximum Entropy Production Principle and Maximum En- 
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tropy Principle still is not clear (see [U [2] for an attempt to derive the former 
from the latter and |15| for a comment). This study expresses the view that 
regardless of the veracity of the principle in full generality, it can certainly 
be useful in application in the field of climate sciences, as well as in many 
others. 

Acknowledgements 

The authors would like to thank Jean-Louis Dufresne for fruitful discussions 
and advice, as well as Gilles Ramstein and Didier Roche for their stimulating 
remarks and constant encouragements. 

References 

[1] R. Dewar. Information theory explanation of the fluctuation theo- 
rem, maximum entropy production and self-organized criticality in non- 
equilibrium stationary states. Journal of Physics A: Mathematical and 
General, 36:631-641, 2003. 

[2] R. Dewar. Maximum entropy production and the fluctuation theorem. 
Journal of Physics A: Mathematical and General, 38:371-381, 2005. 

[3] J.-L. Dufresne, R. Fournier, C. Hourdin, and F. Hourdin. Net exchange 
reformulation of radiative transfer in the CO2 15/im band on Mars. J. 
Atmos. Set, 62:3303-3319, 2005. 

[4] W. M. Elsasser. Heat Transfer by Infrared Radiation in the Atmosphere., 
volume 6 of Harvard Meteorological Studies. Harvard University Press, 
1942. 

[5] C. Essex. Radiation and the irreversible thermodynamics of climate. J. 
Atmos. Sci, 41(12):1985-1991, 1984. 

[6] C. Essex. Radiation and the violation of bilinearity in the thermodynam- 
ics of irreversible processes. Planetary and Space Science, 32(8):1035- 
1043, 1984. 

[7] V. Eymet, J.-L. Dufresne, P. Ricchiazzi, R. Fournier, and S. Blanco. 
Long-wave radiative analysis of cloudy scattering atmospheres using a 
net exchange formulation. Atmospheric Research, 72:239-261, 2004. 



16 



[8] T. Fichefet and M. Maqueda. Sensitivity of a global sea ice model to the 
treatment of ice thermodynamics and dynamics. Journal of Geophysical 
Research, 102:12609, 1997. 

[9] T. Fichefet and M. Maqueda. Modelling the influence of snow accu- 
mulation and snow-ice formation on the seasonal cycle of the antarctic 
sea-ice cover. Climate Dynamics, 15:251-268, 1999. 

[10] K. Fraedrich and F. Lunkeit. Diagnosing the entropy budget of a climate 
model. Tellus A, 60(5):921-931, 2008. 

[11] J. Gerard, D. Delcourt, and L. Francois. The maximum entropy pro- 
duction principle in climate models: application to the faint young sun 
paradox. Q. J. R. Meteorol. Soc, 116:1123-1132, 1990. 

[12] R. Goody. A statistical model for water vapour absorption. Q. J. R. 
Meteorol. Soc, 78:165-169, 1952. 

[13] H. Grassl. The climate at maximum entropy production by meridional 
atmospheric and oceanic heat fluxes. Q. J. R. Meteorol. Soc, 107:153- 
166, 1981. 

[14] J. Green. Division of radiative streams into internal transfer and coohng 
to space. Q. J. R. Meteorol. Soc, 93:371-372, 1967. 

[15] G. Grinstein and R. Linsker. Comments on a derivation and application 
of the 'maximum entropy production' principle. Journal of Physics A: 
Mathematical and Theoretical, 40:9717-9720, 2007. 

[16] F. Hourdin, I. Musat, S. Bony, P. Braconnot, F. Codron, J.-L. Dufresne, 
L. Fairhead, M.-A. Filiberti, P. Friedlingstein, J.-Y. Grandpeix, G. Krin- 
ner, P. LeVan, Z.-X. Li, and F. Lott. The LMDZ4 general circulation 
model: climate performance and sensitivity to parametrized physics with 
emphasis on tropical convection. Climate Dynamics, 27:787-813, 2006. 

[17] E. Jaynes. Information theory and statistical mechanics. Physical Re- 
view, 106:620-630, 1957. 

[18] J. Kiehl and K. Trenberth. Earth's annual global mean energy budget. 
Bull. Amer. Meteor. Soc, 78:197-208, 1997. 

[19] A. Kleidon, K. Fraedrich, T. Kunz, and F. Lunkeit. The atmospheric 
circulation and states of maximum entropy production. Geophys. Res. 
Lett, 30(23) :2223, 2003. 



17 



[20] A. Kleidon and R. Lorenz, editors. Non- equilibrium Thermodynamics 
and the Production of Entropy: Life, Earth, and Beyond. Springer, 
Berlin, 2005. 

[21] A. Lacis and J. Hansen. A parameterization for the absorption of solar 
radiation in the Earth's atmosphere. J. Atmos. Sci, 31:118-133, 1974. 

[22] J. Li. Gaussian quadrature and its application to infrared radiation. J. 
Atmos. Sci, 57:753-765, 2000. 

[23] C. Lin. An extremal principle for a one-dimensional climate model. 
Geophys. Res. Lett, 9:716-718, 1982. 

[24] E. N. Lorenz. Generation of available potential energy and the intensity 
of the general circulation. Dynamics of Climate, pages 86-92, 1960. 

[25] R. Lorenz, J. Lunine, P. Withers, and C. McKay. Titan, Mars and Earth: 
Entropy production by latitudinal heat transport. Geophys. Res. Lett, 
28(3):415-418, 2001. 

[26] V. Lucarini. Thermodynamic efficiency and entropy production in the 
climate system. Phys. Rev. E, 80:021118, 2009. 

[27] V. Lucarini, K. Fraedrich, and F. Lunkeit. Thermodynamic analysis 
of snowball Earth hysteresis experiment: Efficiency, entropy production 
and irreversibility. Q. J. R. Meteorol. Soc, 136:2-11, 2010. 

[28] V. Lucarini and F. Ragone. Energetics of PCMDI/GMIP3 climate mod- 
els: Net energy balance and meridional enthalpy transports. Rev. Geo- 
phys., 49:RG1001, 2011. 

[29] G. Madec, P. Delecluse, M. Imbard, and C. Levy. CPA version 8.1 
Ocean general circulation model reference manual, 3. Technical report, 
LODYC, 1997. 

[30] O. Marti, P. Braconnot, J.-L. Dufresne, J. Bcllicr, R. Benshila, 
S. Bony, P. Brockmann, P. Cadule, A. Caubel, F. Codron, N. D. 
Noblet-Decoudre, S. Denvil, L. Fairhead, T. Fichefet, M.-A. Foujols, 
P. Friedhngstein, H. Goosse, J.-Y. Grandpeix, E. Guilyardi, F. Hourdin, 
A. Idclkadi, M. Kageyama, G. Krinner, C. Levy, G. Madec, J. Mignot, 
1. Musat, D. Swingedouw, and C. Talandier. Key features of the IPSL 
ocean atmosphere model and its sensitivity to atmospheric resolution. 
Climate Dynamics, 34:1-26, 2010. 



18 



L. Martyushev and V. Seleznev. Maximum entropy production principle 
in physics, chemistry and biology. Physics Reports, 426:1-45, 2006. 

R. McClatchey, J. Selby, F. Volz, R. Fenn, and J. Garing. Optical prop- 
erties of the atmosphere. Air Force Cambridge Res. Lab., 1972. 

D. O'Brien and G. Stephens. Entropy and climate, ii: Simple models. 
Q. J. R. Meteorol. Soc, 121:1773-1796, 1995. 

H. Ozawa and A. Ohmura. Thermodynamics of a global-mean state 
of the atmosphere — a state of maximum entropy increase. J. Climate, 
10:441-445, 1997. 

G. Paltridge. Global dynamics and climate-a system of minimum en- 
tropy exchange. Q. J. R. Meteorol. Soc, 101:475-484, 1975. 

G. Paltridge. The steady-state format of global climate. Q. J. R. Me- 
teorol. Soc, 104:927-945, 1978. 

S. Pascale, J. Gregory, M. Ambaum, and R. Taillcux. Climate entropy 
budget of the hadcmS atmosphere-ocean general circulation model and 
of famous, its low-resolution version. Climate Dynamics, pages 1-18, 
2009. 

O. Pauluis and 1. Held. Entropy budget of an atmosphere in radiative- 
convective equilibrium, part i: Maximum work and frictional dissipation. 
J. Atmos. Sci, 59(2):125-139, 2002. 

O. Pauluis and I. Held. Entropy budget of an atmosphere in radiative- 
convective equilibrium, part ii: Latent heat transport and moist pro- 
cesses. J. Atmos. Sci, 59 (2): 140-149, 2002. 

J. P. Peixoto and A. H. Oort. Physics of Climate. Springer: New- York, 
1992. 

J. P. Peixoto, A. H. Oort, M. D. Almeida, and A. Tome. Entropy budget 
of the atmosphere. J. Ceophys. Res, 96:10981-10988, 1991. 

W. Peltier. Global glacial isostasy and the surface of the ice-age Earth: 
the ICE-5G (VM2) model and GRACE. Annual Reviews of Earth and 
Planetary Science, 32:111-149, 2004. 

T. Pujol and J. Fort. States of maximum entropy production in a one- 
dimensional vertical model with convective adjustment. Tellus, 54:363- 
369, 2002. 



19 



[44] T. Pujol and J. Llebot. Extremal climatic states simulated by a 2- 
dimensional model: Part i: Sensitivity of the model and present state. 
Tellus A, 52:422-439, 2000. 

[45] C. Rodgers and C. Walshaw. The computation of infra-red cooling rate 
in planetary atmospheres. Q. J. R. Meteorol. Soc, 92:57-92, 1966. 

[46] G. Stephens. The parameterization of radiation for numerical weather 
prediction and climate models. Mon. Wea. Rev, 112:826-867, 1984. 

[47] G. Stephens and D. O'Brien. Entropy and climate, i: Erbe observations 
of the entropy production of the earth. Q. J. R. Meteorol. Soc, 119:121- 
152, 1993. 

[48] K. Trenberth and J. Caron. Estimates of meridional atmosphere and 
ocean heat transports. J. Climate, 14:3433-3443, 2001. 

[49] K. Trenberth, J. Fasullo, and J. Kichl. Earth's global energy budget. 
Bull. Amer. Meteor. Soc, 90(3):311-323, 2009. 

[50] K. Trenberth and A. Solomon. The global heat balance: Heat transports 
in the atmosphere and ocean. Climate Dynamics, 10:107-134, 1994. 

[51] S. Valcke. 0ASIS3 user guide (prism_2-5). Technical report, PRISM 
support initiative, 2006. 

[52] A. Vos and P. Wei. The efficiency of the conversion of solar energy 
into wind energy by means of Hadley cells. Theoretical and Applied 
Climatology, 46:193-202, 1993. 

[53] W. Wu and Y. Liu. Radiation entropy flux and entropy production of 
the earth system. Rev. Geophys., 48:RG2003, 2010. 

[54] P. Wyant, A. Mongroo, and S. Hameed. Determination of the heat- 
transport coefficient in energy-balance climate models by cxtremization 
of entropy production. J. Atmos. Sci, 45(2):189-193, 1988. 



20 



Profile UH,o {g.cm-') T (K) T^^^^ (K) 



Sub- Artie 
Mid-Latitude 
Tropieal 



272 0.432 295.6 
283 0.436 295.9 
300. 0.444 296.5 



1.5 
2.3 
5.1 



Table 1: Column-integrated water vapour density {uh2o)i characteristic tem- 
perature (T) and resulting longwave atmospheric emissivity [to) for the three 
standard atmospheric profiles. The last column {T^^^^) gives the global 
mean temperature in the MEP state for pre-industrial surface albedo, with 
uniform atmospheric emissivity. 



Peixoto et al PLASIM HadCM3 MEP Pre-industrial MEP LGM 
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17 
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atotai 892 883 912 906 892 



Table 2: Comparison of the entropy budget of the Earth computed from 
observations ([IT], some values have been updated after [20j) and two models 
of different complexity (PLASIM: [10], HadCM3: [STJ), in mW.m-'^.K-\ 
ares designates all the other contributions to the total entropy production 
rate, including eventual imbalance. 
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Figure 1: Energy exchanges in one model box : incoming solar energy (thin 
solid lines) absorbed by atmosphere '^aY ^^"^ ground \&^]^, longwave radia- 
tive exchanges (dashed) ; exchange between the surface and the atmosphere 
'^^gai cooling to space of the surface and of the atmosphere '^l^) and 
large scale heat transport convergence {Q and surface sensible heat (g) (thick 
solid arrows). 
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Figure 2: Surface albedo forcing for the MEP model, computed from the 
output of the IPSL_CM4 simulations. Left : annual mean surface albedo 
values from the pre-industrial simulation with the IPSL model. Right : an- 
nual mean surface albedo values from the LGM ice extent simulation with 
the IPSL model. The surface albedo is higher over continents than over the 
open ocean, and it is higher in high-latitude areas than in low-latitude areas 
because of the presence of snow and ice. The main difference between pre- 
industrial and LGM surface albedo distributions is the presence of a large 
ice-sheet over North America and Scandinavia during the LGM. 




Figure 3: Results for surface temperature. Left: Surface temperature at MEP 
state for present conditions, in ° C. Right: Surface temperature difference 
between the MEP state and the IPSL simulation. Contour lines interval is 
10 ° C, positive contours are drawn in solid lines, negative contours in dashed 
lines and the null contour clS Qb dotted line. 



23 




Figure 4: Left: Total Meridional energy transport predicted by MEP for 
present conditions (solid) and for the IPSL model (dashed), in PW. Right: 
Net radiative budget at the top of atmosphere as a function of latitude (zon- 
ally averaged) for the MEP state (solid) and for the IPSL model (dashed) at 
present conditions, in W.m~'^. 




Figure 5: Left: Surface heat flux predicted by MEP for present conditions, 
in W.m~'^. The mean value is 55W.m~'^. Right: Difference between the 
surface energy flux at MEP state and the surface sensible heat flux for the 
IPSL simulation for present conditions. Contour lines interval is 30W.m~'^, 
positive contours are drawn in solid lines, negative contours in dashed lines 
and the null contour as a dotted line. 
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Figure 6: Pre-industrial surface temperature difference between the MEP 
model with clouds, including a topography correction, and the IPSL model. 
Contour lines represent the -5 isoline (dashed), the null isoline (dotted) and 
the +5 isoline (solid). This figure is to be compared with figure [sl right. 




Figure 7: Left: Surface temperature difference between Last Glacial Maxi- 
mum conditions and present conditions at MEP states. Right: Surface tem- 
perature difference between Last Glacial Maximum conditions and present 
conditions for the IPSL model. Contour lines space is 10 ° C, positive con- 
tours are drawn in solid lines, negative contours in dashed lines and the null 
contour as a dotted line. The global mean LGM cooling is —2 ° C for the 
MEP model and -2.5 ° C for the IPSL model. 



25 




Latitude Latitude 



Figure 8: Total Meridional energy transport predicted by MEP for LGM 
conditions (solid line) and present conditions (dashed line) at MEP state 
(left) and for the IPSL simulation (right), in PW. 
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